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Abstract 

The transitional boundary layer flow over a flat plate is investigated. The boundary layer flow 
is known to develop unstable Tollmien-Schlichting waves above a critical value of the Reynolds 
number. However, it is also known that this transition can be observed for sub-critical Reynolds 
numbers. In that case the basin of attraction of the laminar state coexists with the sustained 
turbulence. In this article, the trajectory on the separatrix between these two states is simulated. 
The state on the separatrix is independent from the initial condition and is dynamically connected 
to both the laminar flow and the turbulence. Such an edge state provides information regarding the 
basic features of the transitional flow. The solution takes the form of a low speed streak, flanked 
by two quasi-streamwise sinuous vortices. The shape of the streaks is close to that simulated with 
the linear optimal perturbation method. This solution is compared to existing results concerning 
streak breakdown. The simulations are realized in a temporal framework for a local boundary 
layer, with periodic boundary conditions in the streamwise direction. A dedicated model, based 
on a scale separation, is presented. The mean flow is a solution of the Prandtl boundary layer 
equations while the superposed small-scale fluctuations are a solution of the periodic Navier-Stokes 
equations. The model is validated with turbulent flow simulations and satisfactorily reproduces 
the physical characteristics of a boundary layer flow, especially in the outer region, where external 
fluid is entrained toward the boundary layer. 



Electronic address: damien.biau@univ-poitiers.fr 



I. INTRODUCTION 



The boundary layer transition to turbulence is a long standing subject of research since 
it presents a good compromise between academic and industrial configurations. The delay 
of transition may be of applied interest, for example, to overcome turbulent aerodynamic 
drag. On the other hand, promoting transition permits an increase in mixing, heat transfer 
and other exchange processes. Since the first results by Lord Rayleigh, more than a century 
ago, this theoretical subject remains fascinating in its richness, as attested to by the regular 
number of related publications. It can be regarded as a benchmark problem for stability 
analysis. This simple shear fiow also raises issues associated with, for instance, receptivity, 
chaos in open fiows, modal versus algebraic growth, the emergence and the sustaining of 
coherent structures. 

The control parameter governing boundary layer transition is the Reynolds number, de- 
noted Re, which can be defined in different ways. Hereafter the characteristic quantities 
are the boundary-layer displacement thickness 6* and the free-stream velocity Uoo, then 
Res* = Uao^* /v. The minimal Reynolds number below which disturbances decrease mono- 
tonically is found to be Res* = 17 4^. The laminar state is globally stable, any finite- 
amplitude perturbations are monotonically decreasing. For Res* > 519.4, the laminar state 
is no longer an attractor and Tollmien-Schlichting (hereafter denoted to as T-S) travelling 
waves are amplified. In the range 17 < Res* < 519.4, the fiow is conditionally stable, its 
stability depends on the disturbance shape and energy. Infinitesimal perturbations can be 
transiently amplified and some coherent states, or traveling wave solutions, are expected 
1| reported the possibility of sustained turbulence up to the minimum 



to cohabit. Spalart 
Reynolds number Res* = 400 



491 . This is characteristic of the the sub-critical nature of the 



instability, indeed the nonlinear Tollmien-Schlichting waves, i.e. the saturated T-S waves 
with a finite amplitude, can be prolonged below the critical value. For these two-dimensional 
equilibrium states Koch J2!] found the nonlinear critical Reynolds number: Res* = 500. 

In addition to these non-linear aspects for the sub-critical transition, an important step 
was made, three decades ago, with the discovery of a possible transient growth for the 
perturbation energy in a linear framework. While the classical analysis focuses on the the 
asymptotic behavior determined by the least stable mode, a general, small disturbance is 
in fact a weighted combination of linear eigenfunctions. Because of the non-normality of 
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the linearized stability operator, there is the potential for very large transient amplification 
of the disturbance energy, even in nominally stable flow conditions. For unstable flows the 
transient growth is able to overcome the classical T-S instability, in this case the transition 
is called bypass. This scenario is encountered for boundary layers subjected to moderate 
perturbations, for example a free stream turbulence with levels comprised between 1 and 
6% of the external velocity. Unlike the single mode case, the growth property depends on 
the disturbance environment, like surface roughness or free-stream turbulence. To overcome 
this difficulty the initial receptivity process can be studied using the optimal perturbation 
method : the inflow or initial condition is obtained through the maximisation of the transient 
energy growth. For the boundary layer, the temporal analysis was flrst performed by Butler 
& Farrell |3[ and extended later to the spatial non-parallel boundary layer by Andersson et 
al. [4 1 and Luchini These works have shown the same features, the optimal perturbations 
take the form of damped streamwise vortices which generate streamwise streaks through a 
lift-up effect. The ratio between the input energy (vortices) and the output (streaks) is 
proportional to the Reynolds number: it is for this reason that the process is relevant to 
the bypass transition. Another bypass scenario is triggered by a pair of oblique waves, 
see Schmid & Henningson [6] and Berlin et al. [7]. In this case, the streamwise vortices 
are generated by quadratic interaction between the oblique wave and, in turn, form the 
streamwise streaks. 

The linear mechanisms described so far are not sufficient to trigger transition in the 
sub-critical regime. In all cases, if the perturbations do not exceed some threshold, they 
eventually decrease and the flow state stays in the laminar basin of attraction. Thus the 
linear analysis has to be pursued by a weakly nonlinear or secondary stability analysis. 
The secondary instabilities can be differentiated following their spanwise symmetry, distin- 
guishing the varicose (symmetric) from the sinuous (antisymmetric) structures. Again the 
streaks and T-S waves differ in their development. The nonlinear development of T-S waves 
instigates three dimensionality, giving rise to longitudinal streaks and varicose A-shaped 
vortices. These vortices can be aligned in the streamwise direction (K-type as reference to 
Klebanoff et al. [8]) or staggered (H-type as reference to Herbert ^9)]). The late stage of 
oblique waves transition exhibits similar A structures, see Berlin et al. [t]. More details 
can be found in the book by Schmid & Henningson [lO|. On the other hand, during the 
streak breakdown, computed by direct numerical simulations by Brandt et al. both 
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sinuous and varicose modes are observed, with the former being more hkely to occur. This 
feature is in agreement with various analyses: a modal inviscid stability by Andersson et 
al. non-modal growth by Hoepffner et al. [l^ and experimentally by Swearingen & 

Blackwelder 1^. While these studies show similar conclusions, it should be noted that the 
streaks considered are generated in different ways, by free-stream turbulence, as a result of 
optimal perturbations or through a Gortler instability. For further details, a recent overview 
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of streak breakdown in bypass transition can be found in Schlatter et al. 

A commonly observed feature at the onset of wall turbulence is the development of 
streamwise low-speed streaks. More precisely, the generic transitional state is composed of 
three main elements: the streamwise rolls sustain the streamwise streaks, by redistributing 
the mean momentum of the mean shear flow. Then the streaks become the support for 
wake-like instabilities, due to the spanwise or wall-normal shear, via modal or transient 
growth mechanisms, see Schoppa & Hussain ij]. Finally the nonlinear (quadratic) self- 
interaction of that streamwise fluctuations directly regenerates the streamwise rolls, closing 
the wall-cycle loop. 

For developed turbulence, such constituents are involved in a self-sustaining process and 
they represent the minimal self-sustaining flow unit capable of surviving viscous decay, as 
demonstrated for the plane Poiseuille flow by Jimenez & Moin 17|. The minimal flow unit 
corresponds to the smallest computational box in which the turbulence may be sustained. 
For each wall layer, a single set of coherent structures was identified, composed of a pair 
of high and low-speed streaks with a pair of counter-rotating, quasi-streamwise vortices. 
The application to the plane Couette flow allowed Hamilton et al. [18] to define the self- 
sustaining process as a quasi-cyclic dynamic, comprising quiescent periods interspersed by 
bursts of shorter duration. 

While the coherent structures were previously obtained by continuation of a sub-critical 
instability by Koch [3| , or by continuation from an unstable brother flow (for example from 
Taylor- Couette to plane Couette flow, see Nagata [19|), the paradigm of the self-sustaining 
process provides a method by which to calculate the exact coherent structures by contin- 
uation of solutions that bifurcate from a streaky flow as explained in Waleffe [20]. Other 
such solutions were later found for plane channel flow and pipe flows. As for the solution 
continued from the marginal instability, they survive down to Reynolds numbers below the 
critical value for turbulence onset. These results suggest that the underlying process is 
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generic and fundamental to both developed turbulence and transition, see Wang et al. \21 \. 



The contribution of Waleffe 



22l | should be mentioned here: he described the unity between 



turbulence in the minimal flow unit, the self-sustaining process of the developed turbulence 
and the role of the embedded traveling wave solutions. It is not the aim here to provide 
a complete bibliography; overviews are provided by Cvitanovic & Gibson [23], Kawahara 



2J] and Panton 25|]. Beyond the various theoretical approaches, a promising and exciting 



result i_s_the direct observation of these traveling waves in the pipe flow experiments by Hof 
et al. 



26l |. For spatially developing flows, Duriez et al. [27| demonstrated, experimentally, 



the relevance of the self-sustaining process in a boundary layer. By varying the amplitude 
of a vortex generator, i.e. small cylinders, the results show a threshold above which the 
streamwise vortices are regenerated from the destabilization of the streaks. 

The minimal flow unit was introduced by Jimenez & Moin [171 ^ ^'^'^^ ^'^^ studying 
near wall turbulent flow with a lower complexity. In the sub-critical case, another strategy 
for minimizing chaotic activity consists in maintaining the trajectory on the separatrix, 
namely the surface separating the laminar basin of attraction and the turbulent trajectories. 
A set of initial states that dynamically evolve to the laminar attractor forms its basin of 
attraction, which is separated from the rest of the state space by the separatrix. With a 
repeated bisection method it is possible to maintain the transitional state, for arbitrary long 
times, on the separatrix. This method was initially proposed by Toh & Itano in 2003 28| . 
who computed a periodic-like trajectory of channel flow. The edge state maintained on the 
separatrix is called the edge of chaos by Skufca et al. [29|. The method was applied latter 
by Schneider et al. [30] to a pipe flow and by Schneider et al. [31[ to a plane Couette flow. 
It is expected that knowledge of these finite-amplitude solutions will contribute much to the 
unraveling of the structure of transitional turbulent flows. 

In this article we are interested in the transitional boundary layer for the sub-critical 
regime. The main objective is to compute a sustained solution, dynamically connected to 
the laminar and the turbulent states. The flow configuration is chosen such that the two 
attractors coexist and the trajectory is maintained on the separatrix with a bisection method. 
The results are compared to the various linear and nonlinear scenarios, mentioned above, 
describing the sub-critical laminar-turbulent transition. Whereas the path to turbulence 
depends on the shape and the amplitude of the initial condition, in contrast, the converged 
trajectory on the separatrix does not depend on its initial condition. Thus the corresponding 
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edge state can provide some general features of transitional states. 

A secondary objective concerns the numerical issues associated with the simulation of 
boundary layer flows, since the previous work for simulating the separatrix were realized for 
parallel flows. In order to extend the method to the nearly parallel boundary layer flow, it is 
necessary to keep the periodic condition in the streamwise direction. However, a temporal 
boundary layer is not satisfactory since the boundary layer thickness increases in time. In 
the next section a model is presented to simulate a local boundary layer, namely with a 
constant displacement thickness (5*); this model owes much to the previous work by Spalart 



II. METHOD 



The overall goal of this section is to obtain a model, with periodic boundary condition in 
the streamwise direction, which can provide an approximate solution of the local boundary 
layer. The coordinates and velocities are made non-dimensional, respectively, using the 
boundary- layer displacement thickness 6* = 1 — U dy and the free-stream velocity Uoq. 
The associated Reynolds number is Res* = Uoo5* /v. 

To retain the advantages of the spectral Fourier decomposition in the streamwise direc- 
tion, it is possible to add a fringe region downstream of the physical domain. The outflowing 
fluid is constrained, via a volume force (which is signiflcantly non-zero only within the fringe 
region), to the desired inflow condition, see for example Brandt et al. |32|. The problem 
of specifying the inflow boundary condition remains. Lund et al. suggested that the 
velocity fleld at a downstream station be rescaled and re-introduced at the inlet. The rescal- 
ing is based on the similarity laws of the boundary layer, the difficulty comes from the two 
different wall-normal length scales : the law of the wall in the inner part and the defect 
law in the outer part of the boundary layer. This two-layer nature of the turbulence intro- 
duces spurious periodicity into the time series during the recycling (see Siemens et al. |34|). 
As a consequence, a substantial initial part of the computational box has to be discarded. 
The numerical cost leads us to reject this method in our case. An alternative consists in 
simulating a temporal boundary layer in a moving frame of reference, see Spalart & Yang 



351]. However the undisturbed boundary layer is not a solution of the unsteady boundary 



layer equations. In addition, the boundary layer thickness increases in time so the base flow 
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eventually becomes unstable, thus it is not possible to reconcile the existence of an edge 
state with the large time necessary to obtain converged results. Thus it appears necessary 
to introduce a specific method for a local turbulent boundary layer. Following the method 
introduced by Spalart the velocity field u is decomposed as a local boundary layer profile 
U plus the deviation: 

u(x, y, z, t) = U(x, y, t) + u'{x, y, z, t). (1) 

The mean fiow is solution of the turbulent boundary layer equations: 
dU 9V _^ 
dx dy 

dU ( dU dU\ , 1 d'^U dvh? ^ ' 

I Re 



(9x dy) 



dt \ dx dy J Re dy"^ dy 

The overbar on the Reynolds stress indicates a streamwise and spanwise averaging. The 
momentum equation is Reynolds number dependent because the time scale is of order 6* /Uoo, 
in agreement with the time scale of the streamwise fiuctuations (see equations [6]). One major 
consequence of the non parallel effect is the continuous entry of external irrotational flow 
while the boundary layer thickness increases downstream. To take into account this effect, 
a new coordinate is introduced : t] = y/6* and the contravariant velocities (letters topped 
by a tilde) are used: 

U = U 

dU dU dr] dU 

— = \ — - y — (3) 

dx dx dx dy 

dx 

The wall-normal derivatives, dU /dy and dV /dy are computed directly from U and V . This 
transformation is classical for adapted meshes, see also Spalart . With this transformation 
V is zero in the wall layer and becomes negative farther from the wall, thus the non-turbulent, 
or irrotational, fluid is entrained into the boundary layer. 

The boundary layer equations are solved between two locations: Xi and X2, for an 
inflnitesimal distance X2 — Xi. In this sense the spatial boundary layer is simulated locally, 
for a given Re^*. This distinguishes the present model from the temporal boundary layer 
flow. The former describes, locally, the temporal evolution of a spatial boundary layer [i.e. 
fixed 5*) while the latter describes a temporal boundary layer fiow {i.e. with 5* evolving in 
time) . 



At the location Xi, the displacement thickness S*{Xi) is determined, through a shooting 
method, such that 6*{X2) = 1. Now that r] is defined, the streamwise velocity at Xi is 
determined, assuming a self-similar behavior : 

U{X,, T]) = U{X2, y/5*{X,)) (4) 

It is now possible to compute the streamwise derivatives, in equations [2] and [3l by taking 
the differences : 

df] _ y-y/6*{X,) 

dx X2 — Xi /_s 
dU _ U{X,, y/6*{X,))-U{X2, y) ^ ' 

dx X2 - Xi 

The results appear to be weakly sensitive to the spacing X2 — Xi provided it is sufficiently 

small, at least lower than 0.5. In the following, the simulations are realized with X2 — Xi = 

10-2. 

The deviation is assumed to evolve on a short scale around a weakly nonparallel mean 
ffow evolving on a long scale. The equation for the deviation truncated at order 1/Res* 
becomes: 



du' ^ dv' ^ dw' ^ 
dx dy dz 

dt \ / dx \ / dy dy dx dz dx 

dv' f~ ,\ dv' frri-r. ,\ dv' dV/Re , ,dv' dp' . 

— + f/ + m' — + (V Re + v') — + v' + w'— = + Av' Re 

dt \ / dx \ ^ dy dy dz dy 

dw' ( ~ A dw' f-Tr/j:, i\ dw' ,dw' dp' 



(~ A aw f-TriT^ ,\ ow ,dw dp . 
+ [U + u') —+ (V Re + v') — + w'— = -^ + Aw' Re 
V / dx \ ^ dy dz dz 

(6) 

with a compact form of the Laplacian operator A = d'^ jdx^ + d'^ jdy^ + d"^ /dz^. These 
equations are associated with homogeneous boundary conditions at the wall and at the upper 
boundary. The results appear not to be inffuenced by the external homogeneous boundary 
condition since the simulated ffuctuations are irrotational and behave asymptotically as 



exp(— A; y) with k = a/o^T^. 

The two strong assumptions are the periodicity applied to streamwise {i.e. x-independent) 
streaks and vortices in equation [6] and the rescaling of the mean streamwise velocity. On the 
other hand, the model is self-governing, which means there is no need for an empirical ad- 
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justment with previously existing data. The model could be interesting in order to generate 
an inflow condition for a spatially developing turbulent boundary layer. 

The equations are discretized using a Chebyshev collocation method in the wall-normal 
direction {y) and Fourier in the streamwise (x) and spanwise {z) direction. In the wall- 
normal direction, the Gauss-Lobatto-Chebyshev points {() are expanded from [—1, 1] to 
[0, yoo] throughout an algebraic mapping: y = (?/oo^)/2 (1 + C)/(l + ^C)- Where h = 0.5 
is a stretching factor and yoo = 15 is the location of the upper boundary. The differenti- 
ation operators are suitably adjusted by the analytical Jacobian of the mapping. Integral 
quantities are obtained using similarly-adjusted Gaussian quadrature rules. For the time- 
integration an order three semi-implicit finite difference scheme is used, associated with 
a prediction/projection method. In order to compute a pressure unpolluted by spurious 
modes, the pressure is approximated by polynomials {Pn~2) of two units lower-order than 
for the velocity (-P/v), see Botella et al. 



The simulation must resolve the entire range of scales of the turbulent motion : from the 
boundary layer thickness down to the dissipative scale. As an indication, the mean spanwise 
spacing of the streaks in the buffer layer {y^ < 30) is about 100 wall units, the superscript + 
designates quantities expressed in the wall units based on the kinematic viscosity v and the 
skin friction velocity Ur- The first constraint implies a large domain while the latter imposes 
a large number of grid points. The period in the streamwise and spanwise directions are 
Lx = 50 and Lz = 15. In the transverse directions the grid spacing is Az^ ^ 6 and 
Ax+ ~ 20, while in the wall-normal direction, at least 14 points are within 9 wall units of 

n 

the wall. These requirements match those by Spalart pj], except that the domain is slightly 
shorter in the present simulations. The associated number of points in the x-, y-, and z- 
directions is (64 x 75 x 64) at Res* = 500 and (128 x 95 x 128) at Res* = 1000, while the 
timesteps are At = 0.1 and At = 0.05, respectively. 

For zero fluctuations, the laminar Blasius flow is reproduced with an absolute error lower 
than 0.15%. The model is also able to correctly reproduce the linear growth of ToUmien- 
Schlichting waves. The ability of the model to reproduce the turbulent boundary layer now 
needs to be demonstrated by comparison with other results. The mean velocity profiles are 
displayed in figure [H 

The agreement is correct in the buffer layer (?/"•" < 30), an area in which the DNS is very 
reliable. Just above, in the log layer, the agreement is also satisfactory but farther from 
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FIG. 1: Mean velocity profiles at Res* = 500 (bottom curves) and Res* = 1000 (upper curves). 
Continuous lines : present results, dashed lines Spalart [l|]. The dotted line indicate the classical 
viscous sublayer and the log law k^^ lny+ + C, with k = 0.4 and C = 5. 

the wall, in the wake region, the flow is misrepresented as expected from a nearly parallel 
model. The skin friction velocity is overestimated for both nearly parallel simulations. 




400 



FIG. 2: Reynolds stress distribution, normalized by the friction velocity, at Res* = 1000. Contin- 
uous line : present results, dashed line Spalart [l|, dotted line Schlatter &i Orlii 
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This effect is more pronounced for the Reynolds shear stress (r = —uv~^), displayed 
in figure |2l which is wider than expected. The slight differences between present results 
and those by Spalart jlj] can be imputed to the stronger parallel-flow hypothesis since the 
streamwise streaks and rolls are rescaled in the results by Spalart l|. 

To conclude with quantitative comparisons, the Reynolds number based on the momen- 
tum boundary layer thickness is Ree = 307, for the case Res* = 500, which is in agreement 
with the value Ree = 303 computed by Spalart it- For the case Res* = 1000 the present 
result is Ree = 683, in relatively good agreement with the values 677 found by Spalart [l| 
and 670 found by Schlatter et al. [37]. The fluctuating wall-shear stress are r^^rms = 0.372 
and 0.3922 for Res* = 500 and 1000, respectively. These values are close to the classical 
value of 0.4, also in good agreement with the relation given by Alfredsson et al. jssi and 
Orlii & Schlatter Isoj: T^,rms = 0.298 + 0.018 lni?e^ leading to: r^,^^, = 0.388 and 0.401. 

However, we will see that the separatrix is close to the laminar state, in particular the 
Reynolds shear stress for the edge state is three time smaller than for the turbulent state. 
In addition, the turbulent inner layer being absent, the mean flow U scales exactly with the 
outer scale 6*, in agreement with the model (see equation H]) . Bearing in mind that the main 
issue of the paper is the understanding of the edge state and not the prediction of turbulent 
state, these results are deemed sufficient validation of the code for the present purpose. 

A trajectory on the separatrix is bounded by initial conditions which eventually decay 
or become turbulent. Practically, we consider a turbulent velocity field u, decomposed as 
shown in equation [H with weighted fluctuations: 

u(a;, y, z, t) = \J{x, y, t) + Au'(a;, y, z, t). (7) 

For values of the weight A close to unity, the flow stays turbulent. By reducing A, the 
initial condition moves closer to the laminar profile and will not become turbulent again. 
Therefore an interval of A values can be found, bounded at one end by a trajectory that 
becomes turbulent and at the other by a trajectory that returns to the laminar state. In 
between two such trajectories lies one which resides for a substantial time interval on the 
separatrix. Thus the numerical method is based on a bisection algorithm in order to find the 
amplitude of two initial conditions on either side of the edge state. The bisection algorithm 
is initialised with a turbulent flow but the choice of initial state is not critical since arbitrary 
initial conditions appear to converge to the same solution. This edge state is only a partial 
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attractor, equivalent to a saddle point in a phase space. Since it is necessary to repeat this 
procedure (after every 200 time units in our case) to constrain the solution on the separatrix 
and then to obtain a converged solution. After convergence, the coefficient A is very close 
to unity, with an accuracy of five decimals places, so the edge state can be considered as a 
solution of the nonlinear Navier-Stokes equations. 

It should be noted that this method is different from finding a threshold for transition to 
turbulence for a given initial condition. In the latter case the trajectory is found close to the 
separatrix but it is not exactly on the separatrix. Although there are several paths from the 
laminar state to turbulence, depending on the initial condition, the state on the separatrix is 
independent on the initial condition. In other words, the laminar basin of attraction seems 
to be simply connected. 

Even if the basin of attraction of the laminar state disappears for Res* > 519.4, the neutral 
stability curve, in the figure [3l shows that for streamwise length shorter than = 17.5 the 
Blasius boundary layer is unconditionally stable and the laminar state is always surrounded 
by a basin of attraction, for all Reynolds numbers. 




FIG. 3: Neutral stability curve for the Tollmien-Schlichting perturbation : Reynolds number (Res*) 
versus streamwise length (L^), the unstable domain is the gray area. For streamwise length shorter 
than Lx = 17.5 the Blasius boundary layer is unconditionally stable. The three circles represent 
the parameters used for the simulations. 

Thus it becomes possible to investigate the Reynolds number effect, even for super-critical 
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values, by using a streamwise length equal to 17.5. Eventually three sets of parameters have 
been retained. The main simulation was realized with Res* = 500 with = 30. Two 
complementary simulations have been realized with = 17.5 for two Reynolds numbers: 
Res* = 500 and 1000. In all simulations the numerical parameters are A^^; = 64, Ny = 
75, N, = 64 and At = 0.05. 

In order to check the improvement of the model presented in this section, a preliminary 
simulation has been performed with a parallel Blasius solution as a mean flow. The parallel 
assumption leads to numerical issues, the perturbations extend outside the boundary layer 
and eventually fill all the numerical domain, as for a half-channel flow. 

III. RESULTS 

Firstly in order to illustrate the dynamic on the separatrix, the time evolution of the 
shape factor is depicted on the figure HI The shape factor H = 6*/6, with 6 the momentum 
thickness, is often used to determine the nature of the boundary layer flow. The time 
evolution of the shape factor in figure H] indicates that the separatrix is much closer to the 
laminar state {H = 2.59) as compared to the turbulent state {H = 1.63). 

3 

2.5 
H 

2 



■■"0 2000 4000 6000 8000 10000 
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FIG. 4: Shape factor H versus time for Res* = 500 and = 30. Dotted lines indicate the laminar 
and the (time-averaged) turbulent values. The continuous line corresponds to the separatrix, the 
two dashed lines represents the trajectories on either side of the edge state. 

The bisection method is stopped after 6000 time units. After t = 6000 the trajectory 
is ejected away in two possible directions: toward the laminar attractor or toward the 
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turbulence, plotted with dashed lines in the figure HI The laminar and turbulent trajectories 
remain close for 500 time units, confirming that the numerical solution is very close to the 
separatrix. When the trajectory escapes from the separatrix, along the laminar path, the 
streamwise fiuctuations disappear first and the fiow tends to be formed by quasi-straight 
rolls and streaks, which then decay on a (slow) viscous time scale towards the laminar fixed 
point. It should be noted that the discrepancy between the final value {H = 2.58) and the 
expected Blasius value {H = 2.591) can be explained by the approximated Blasius solution 
computed with equations [2l On the other hand, during the transition to turbulence the 
streamwise streaks and vortices move closer and closer to the wall, prefiguring the developed 
turbulence where the streamwise vortices populate the inner region; however the turbulent 
streaks are not an extension of the transitional streaks. Thus the two-layer nature of the 
wall-turbulence emerges: close to the wall, energy production exceeds dissipation, a part 
of this energy is exported far from the wall where dissipation dominates. In between a 
region develops, comparable to the inertial range of scales in the energy cascade of isotropic 
turbulence, where the production balances the dissipation. But this cannot be considered 
as a logarithmic region for this low Reynolds number. 

The case with higher Reynolds number Res* = 1000 and shorter streamwise length = 
17.5 is plotted in figure [51 The dynamic constrained in a short domain seems to be a periodic 
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FIG. 5: Shape factor H versus time for Lx = 17.5 with Res* = 500 (left) and Res* = 1000 (right). 



limit cycle. The period is found to be almost constant and close to 1200. The simulation with 
a lower Reynolds number {Res* = 500) presents the same oscillations with a period of about 



630. However, as stated by Toh & Itano 



, we cannot conclude whether the periodic-like 



14 



solution is a heteroclinic cycle or an exact periodic solution. On the other hand, for pipe 
flow, Duguet et al. [40| indicate that trajectories on the laminar-turbulent boundary could 
be organized around only a few traveling wave solutions and their heteroclinic connections. 



See also Cvitanovic & Gibson 



23| for the relations between the periodic solutions and 



heteroclinic cycle in plane Couette flows. In the present case, the cyclical behavior seems 
to be a numerical artifact induced by a short streamwise length. In the following we focus 
on the results obtained for the larger box = 30 and a sub-critical Reynolds number 
Res* = 500. 

When numerical convergence is reached, the edge state does not depend on its initial 
condition. In that sense it is interesting to compare this state with the last state of transi- 
tional flow presented in various other investigations. Figure E] shows a cross section of the 
flow fleld at t = 2800. Streamwise vortices are associated with low and high speed streaks 
throughout the lift-up process describing the linear interaction between the rolls and the 
mean shear flow. The non-linear interaction between the rolls and the streaks explains why 
the width of the high speed streaks is larger than that of the low speed streaks: the high 
speed streak, induced by a negative vertical velocity, is moved closer to the wall and spreads 
in the spanwise directions. 




FIG. 6: Representation of the streamwise averaged fluctuations at t = 2800 for Res* = 500 and 
Lx = 30. Tlie cross-stream velocity components (vortices) are depicted with vectors while the 
streaks are identified by the grayed areas. 

The most energetic structure is dominated by the low-speed streak associated with a 
pair of strong counter-rotating vortices near the wall, and above a pair of opposite counter- 
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rotating vortices; these are similar to those found in a circular pipe flow by Schneider et al. 



30l | or in square duct by Biau & Bottaro [41]. This structure, with double-layer vortices is 



a tenuous equilibrium between two antagonistic phenomena: the ejection and sweep events. 
The low speed streak, generated by the two near-wall vortices, is stretched in the spanwise 
direction and compressed in the wall-normal direction. It is the location of large, potentially 
inflectional, shear stress. This tenuous equilibrium can explain the periodic-like behavior 
mentioned before. One cycle consists of two typical intervals: a single-streak development 
and its sudden slip in the spanwise direction. The kinetic energy of the fluctuations is 
dissipated gradually in the former case, while it is recovered quickly in the latter case. This 
cycle is present in each case but it is less regular for the case with larger streamwise length. 

The edge state is dominated by two streamwise vortices with an associated low-speed 
streak; as such, the edge-state bears a resemblance to optimal disturbances. The streamwise 



velocity fluctuation profile {ul^^ = 
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dx dz dt) is presented in figure [71 
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FIG. 7: r.m.s profiles for the streamwise velocity : for the edge state at Res* = 500 (continuous 
line) and Res* = 1000 (dashed line), the optimal streaks (dotted line), and experiments (circles). 
The curves have been normalized with their maximum value. 



The result is compared to the profile from the optimal perturbation theory (see Luchini 
[sl and Andersson et al. [4]) and with experimental data by Coupland 42 1 (T3A- test 
case). The experiment represents profiles of r.m.s. velocity perturbations measured in a 
flat-plate boundary layer at location Res* = 1776, perturbed by grid-generated free stream 
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turbulence. The turbulence level is lower than 1% at the leading edge, thus the turbulence 
is triggered by classical unstable Tollmien-Schlichting waves. That background turbulence 
is the reason why the experimental profile does not tend to zero outside the boundary layer. 
The wall normal mode shape of the streamwise Urms edge state is slightly different from the 
optimal perturbation or the values measured in a transitional boundary layer. In particular, 
the maximum value is attained closer to the wall, this can be explained by the two outer 
vortices, shown in the figure [6l which maintain the structure closer to the wall. The time 
averaged amplitude of Urms is 12%?7oo {Res* = 500) and 9%Uoo {Res* = 1000) for the edge 
state, which is close to the respective values 14%f/oo and 13%?7oo observed in the turbulent 
case. In addition, the rms amplitude of the streamwise fluctuations is at least one order 
of magnitude greater than the values for the crossflow fluctuations, in agreement with the 
lift-up mechanism, indeed Luchini [51] demonstrated that the amplitude ratio between the 
streaks and streamwise vortices scales with the Reynolds number. 

Figure [8] provides a snapshot showing of the three dimensional vortical structure identified 
by the Q-criterion. The Q-criterion locates regions where rotation dominates strain in the 
flow, see Hunt et al. [43]. Letting S and Vt denote the symmetric and antisymmetric 
parts of the velocity gradient tensor Vu, the definition of the Q-factor is given by: Q = 
{QijQij — SijSij)/2. A coherent vortex is defined as a region where Q is positive. 

The main structures consist of quasi-streamwise vortices, of alternating sign, in a stag- 
gered pattern located on the flanks of the low-speed streak located in the middle of the box, 
see also the figure El The observed sinuous structures are in agreement with the streaks 
of stability analysis indicating that the sinuous modes preferentially initiate the breakdown 
process. The overview on streaks-breakdown by Schlatter et al. 15|] also confirm the preva- 
lence of the sinuous structures over the varicose structures for transitional flows. 

However, there exist numerous paths from laminar flow to turbulence, depending on 
the initial conditions, their amplitude or shape. In particular, typical varicose structures, 
taking the form of hairpin or horseshoe-shaped vortices, have been observed in a transitional 

n n 

spatial boundary layer by Wu & Moin [44|] and Cherubini et al. [45]. Nonetheless, the edge 
state on the separatrix is found very close to the laminar state and it is characterized 
by fluctuations with relatively low amplitudes. Thus strong varicose structures appearing 
during the transition to turbulence may be considered as subsequent or secondary structures 
appearing on the path to turbulence, beyond the separatrix. Alternatively those structures 
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FIG. 8: Vortical structures on the separatrix at t = 2800 for Res* = 500 and = 30. 



may be on a different trajectory, initialized outside the laminar basin. 

Wu & Moin [44] present instantaneous flow fields in both transitional and turbulent 
regions which are populated by hairpin vortices. The hairpin-shaped vortices, not just 
quasi-symmetric or one legged hairpins, have been observed clearly in both the later parts 
of the by-pass transition and in the fully developed turbulent region. They also observed a 
peculiar bump in the Reynolds stress curves which is slowly damped downstream but still 
visible in the developed part of the turbulence. Wu & Moin 4J] argued that none of the 
previous simulations were of a genuine spatially developing turbulent boundary layer. From 
a different point of view, it could be argued that the presence of hairpin vortices, in this 
case, is a consequence of the periodic forcing used to trigger the transition. After passage of 
the free-stream turbulence pocket, the streamwise flow is straightened and enters into the 
relatively slower-moving turbulent puff. The passage of this high-speed plug flow past the 
slower fluid, that resides near the wall, leads to a strong spanwise sheet of vorticity at the 
interface, away from the wall. The roll-up of this shear layer, through the Kelvin-Helmholtz 
instability, could be interpreted as the head of the hairpin vortices. This is one of the possible 
explanations for the production of hairpin- vortices presented by [25]. This effect has been 
experimentally investigated, in a different context, by Bandyopadhyay et al. [46] for the 
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turbulent puff in transitional pipe flow. 

n 

Cherubini et al. [45[ have simulated the evolution of a nonlinear optimal impulse (lo- 
calized in time and in streamwise direction) in spatially developing flow. They found the 
threshold amplitude above which the disturbances degenerate into a turbulent puff, suggest- 
ing the existence of relative attractors which are formed by a main hairpin vortex placed at 
the leading edge of the packet, followed by a population of quasi-streamwise vortices. Their 
simulations are realized for a super-critical Reynolds number, thus the solution evolves in 
between the classical path to turbulence (involving the slow increase of a T-S wave-packet) 
and the faster path where the T-S amplification is stimulated by the Orr mechanism (see 
Butler & Farrell [sl). This could explain the presence of hairpin vortices which are the 
classical nonlinear development of the T-S instability (see Schmid & Henningson lo|). 



IV. CONCLUSION 



The knowledge of the edge states is important for understanding the sub-critical tran- 
sition to turbulence, providing a complementary tool, beyond the linear or weakly linear 
approaches. At the boundary between order and chaos, we observe the features of both 
states. There is enough order to extract some information about the mechanisms triggering 
the transition, and enough nonlinearity to provide indications about the dynamics of coher- 
ent structures. The vortical structures show a sinuous symmetry, in agreement with streak 
linear stability analysis. However, observations of transitional and turbulent streaks have re- 
vealed the presence of both sinuous and varicose symmetries at breakdown. One may regret 
the limitation of the model for the large scale structures. Because of the strong hypothesis 
of the model, especially the streamwise periodicity, the results cannot be extended to the 
large scale motions. Thus the results regarding the appearance of sinuous vortices cannot be 
considered as definitive and that open question deserves further investigation. Nonetheless 
the possibility to extend the method to unstable flows, providing that the streamwise length 
is short enough, opens up possibilities for the application to a wide variety of shear flows, 
like the mixing layer for example. 

An open question concerns the possible existence of equilibrium states on the separatrix. 
The present challenge is to describe the flow dynamics as a trajectory in phase space mean- 
dering around typical equilibrium states (time-periodic solutions). This idea can be linked 
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to the concept of active and inactive motion presented by Townsend [47|], where the exact 
coherent solutions provide the active part, while the residual (assumed to be uncorrelated 
to the coherent state) contain the inactive motion. These exact coherent solutions, with 
their stable and unstable manifolds, are typically low- dimensional saddle points which col- 
lectively produce a chaotic repeller that is a random switching from one unstable periodic 
solution to another. The complexity of these dynamics necessitates a more complex iterative 
method than the bisection with a single parameter governing the fluctuation amplitude. A 
perspective would be to pursue the method with two parameters, where the second is used 
to constrain the dynamics of the solution on the separatrix or to compute edge states with 
more than two unstable directions. 
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